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Abstract 

We compute that the growth of the origin occupation-time variance up to time t 
in dimension d = 2 with respect to asymmetric simple exclusion in equilibrium with 
density p = 1/2 is in a certain sense at least ilog(logf) for general rates, and at least 
t(log t) 1 / 2 for rates which are asymmetric only in the direction of one of the axes. These 
estimates are consistent with an important conjecture with respect to the transition 
function and variance of "second-class" particles. 

Abbreviated title: Occupation-time variance in 2D asymmetric exclusion process 
with density 1/2 

AMS (2000) subject classifications: Primary 60K35; secondary 60F05. 

1 Introduction and Results 

It is known that the occupation-time variance at the origin up to time t in asymmetric 
exclusion processes in equilibrium is proportional to t times the expected time a second- 
class particle, beginning at the origin, spends at the origin, that is, t L(l—s/t)p s (0, 0)ds 
where p s (0,j) is the second-class particle transition function. Let us now fix the equi- 
librium density p = 1/2 so that the mean of the second-class particle at time t vanishes. 
Recently, it has been argued, as the variance of a second-class particle at time t start- 
ing initially at the origin-^ j 2 Pt(0,j) in this case-is conjectured to be on the order 
i 4 / 3 in d = 1 and proved (for a closely related resolvent quantity) to be at least 
i 5//4 [5], that the the transition function of the second-class particle decays on order 
t~ 2 / 3 in d = 1 (cf. equation (4.8) |T). In d = 2, the second-class particle variance is 
conjectured as 0(i(logt) 2/ ' 3 ) [Q with a proof (for a resolvent quantity when the process 
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rates are asymmetric only in the direction of one of the axes) ^2]. Perhaps by the 
same sort of reasoning as in [7J, one may claim the second-class transition function 
decays as t~ 1 (logt)~ 1 / 3 in d = 2 (cf. equation (12) pQ). Then, the occupation-time 
variance orders should match second-class particle variance orders in both d = 1 and 
2. We mention also these variance orders have connections to fluctuation orders of the 
current across the origin on which there has been much study (cf. [7], [T2]). 

In this note, we show that the occupation-time variance at the origin diverges in 
a sense in d = 2 when density p = 1/2 at least as t log (log t) for general asymmetric 
rates, and at least as i(logi) 1 / 2 when the asymmetry is only in the direction of one of 
the axes (Proposition ll.Hj l and so is consistent with the above discussion. The methods 
are to link occupation-time variances and certain resolvent H_\ norms, and then to 
use some "free-particle" comparisons of H.T. Yau in the style of Bernardin |2j. 

Model. Informally, the simple exclusion process on Z d is a collection of random 
walks which move with jump rates p(i,j) = p(j — i) independently except in that 
jumps to occupied vertices are suppressed. In this article, we will assume that p is 
finite-range such that its symmetrization (p(-) + p{—-))/2 is irreducible. More formally, 
let S = {0, l} z<i be the configuration space where a configuration n = {rji : i G Z d } 
is a collection of "occupation" coordinates where rji = 1 if i is occupied and r\i = 
otherwise. The exclusion process is a Markov process n(t) evolving on £ with generator 

Lf( v ) = ^jpv-iMi-vMfm-m)- 

Here, rf 3 is the configuration obtained from n by interchanging the values at i and j. 
Let also T t denote the associated semi-group. See jS] for more details. 

It is well-known that there is a family of invariant measures {P„ : < p < 1} 
each of which concentrate on configurations of a fixed density p. These measures take 
form as Bernoulli product measures, that is, P p independently places a particle at each 
vertex with probability p. Let E p denote expectation with respect to P p . Denote also 
by (•, - ) p and || • ||o the innerproduct and norm on L 2 (P p ). 

We also note one can compute, with respect to P p , that the adjoint L* is the 
generator of simple exclusion with reversed jump rates p{— ■). 

Problem and Connection to Second-Class Particles. Consider the centered occupa- 
tion time, say, at the origin up to time t, A p {t) = Jq(?/o(s) — p)ds. The problem is to 
compute the variance of A p (t) under the equilibrium P p . Let of = E p [A 2 p {t)\ denote 
the variance. We compute, using stationarity and basic calculations, that 

of = 2 / [ E p [{ m (s) - p)( V0 (0) - p)]duds 
Jo Jo 

= 2 f ( t -s)E p [(r ]0 (s)-p)( m (u)-p)]d S . 
Jo 

To express the kernel further, consider the "basic coupling" of two systems, the first 
starting under ^ ~ P p (-|r/o = 0) and the second under ^ + 5q, that is with an extra 
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particle at the origin. Let (£(t),R(t)) ~ P denote the coupled process where R(t) 
tracks the discrepancy or "second-class" particle. The joint generator is 

r) = Y,p{j-W-W(C j ,r)-m,r)) 
+ E P(-i%- t (f(e- l ' r , r-i)~ fit, r)) 

i 

+ J>(*)(1 - $r+i)(/(e r ' r+< , r + i) - /(£, r)). 

i 

The first sum refers to jumps not including the discrepancy, while the second and third 
sums correspond to jumps of other particles to the discrepancy position and jumps of 
the discrepancy itself. 
We have then 

£p[fa>M-p)(%(o)-p)] 

= p(l - p)P p (vo(s) = % (0) = 1) - P p ( m (s) = l|r/ o (0) = 0) 
= p(l-p)P[R(s)=0] 

which leads to 

f't f'OO 

\miaf/t = lim 2/9(1-/9) / (1 - s/t)P[R(s) = 0]ds = 2/5(1 -p) / P[R{s) = 0]ds; 

t^oo t^co JO ' JO 

the notation earlier in the introduction reads now Pt(0,j) = P(R(t) = j). 

The second-class particle process R(t), with respect to its own history, is not Markov 
except when the jump rate p is symmetric, in which case, it is a symmetric random 
walk. In general, it is highly dependent on the whole system. However, one can roughly 
think of R(t) as some sort of random walk with mean drift (1 — 2p) Y2i ip(i)- This drift 
vanishes exactly when p is either mean-zero (^2iip{i) = 0) or p = 1/2, and so one 
might think the process is recurrent exactly in this case so that 

lim a\ jt = cxd in d < 2 when p mean — zero or when p = 1/2 

t—>oo 

< oo otherwise. 

This has been established in dimensions d > 3, d = 1, and in d = 2 when p ^ 1/2 j3], 
Hli |£|) HI- Still open it seems is to show the variance is superdiffusive in d = 2 when 
P = 1/2- 

Of key interest is also how fast of /t diverges in d < 2 when p mean-zero or p = 1/2. 
In fact, it has been shown that a\ ~ i 3//2 and tlogt in d = 1 and d = 2 respectively 
when p is mean-zero jl], [§]• When p has a drift (X^pW 7^ 0) and p = 1/2, as 
mentioned earlier, of is conjectured to diverge as i 4//3 and t(logt) 2//3 in d = 1 and 
d = 2 respectively. Indeed, a lower bound on order t 5 / 4- has been shown in |2] in 
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d = 1. The main result of this article (Proposition is to compute in d = 2 when p 
has a general drift and p = 1/2 that 

1 f°° 

liminf- — — — / e- xt P(R(t) = 0)dt > 

A^o log(|logA|) J 

or by integrating the second-class transition function twice 

A 2 f°° 

liminf:; — — / e~ xt afdt > 0. 

A^o log(| log A|) Jo 

When the drift ^ ip{i) is in the direction of one of the axes, the same result holds with 
"log | log A|" replaced by "| log Al 1 / 2 ." Clearly o~ 2 /t diverges regardless, and moreover a 
formal Tauberian analogy would suggest that a 2 is at least on order ilog(logi) in the 
general case and i(logi) 1 / 2 in the more special case. 

We mention that some rough upper bounds in d = 1,2 in the "drift" case when 
p = 1/2 are easy to obtain by a comparison with the symmetrized process, namely 
c 2 < ci(p)i 3//2 in d = 1 and < C2(p)tlogt in d = 2. Although well known, we include 
them for completeness in Proposition 11.21 

Variational Formulas. The method of proof does not work with second-class par- 
ticles, but with certain variational formulas for some resolvent quantities. By a local 
function, we mean a function supported on a finite number of coordinates. 

The generator L can be decomposed into symmetric and anti-symmetric parts, 
L = S + A where S = (L + L*)/2 and A = (L-L*)/2. Since L is Mar kovian, S is a non- 
positive operator. Consider the resolvent operator (A — L)^ 1 : L 2 (P p ) — > L 2 (P p ) well 
defined for A > 0-in particular, (A — L) -1 / = J °° e~ Xs T s fds. Since the symmetrization 
of (A-L)- 1 has inverse (A - L*)(X - S') _1 (A - L) = (A - S) + A*(\ — S)~ 1 A, we have 
the variational formula for / local, 

(/, (A - L)- l f) p = sup {2(f, <p) p - (0, (A - S)<f>) p - (A<j>, (A - L^A^),}. 

<j> local 

Now, as A* = —A and A* (A — S)~ 1 A is a non-positive operator, we have the easy 
bound that (/, (A — L)^ 1 f) p is bounded by its "symmetrization," 

{f,{X-L)- l f) p < sup {2(/,0>„-(&(A-S)0) p } 

(f> local 

= </, (A - £)- % (1.1) 

Upper bounds. Well known upperbounds on o~ 2 follow from two statements which 
we include here for completeness. 

Proposition 1.1 There is a universal constant C\ such that 

°1 < C\t{r}o — p, (t _1 - L)~ l r]o — p) p 
< C x t{no - p, (t' 1 - 5)-% - p) p . 
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Proof. The first line is well-known (with a proof found for instance in Lemma 3.9 
0), and the second is (jTT|) . □ 



Proposition 1.2 In d <2, there exists a constant C2 = C2(d,p,p) where for large t, 
(770 - p, (r 1 - S)' 1 ^ - p) p < 



Ci^ft in d = 1 
C 2 log i in d = 2 

and so 6?/ Provosition M.ll of < CiC^t 3 / 2 in d = 1 and CiC2tlogt in d = 2. 

Proof. This is proved in [I] as follows: Write (770 — p, — S 1 ) -1 ^ — p) p = 
ffe-'/tEpifaW - p)(%(0) - p)]ds = p[l - p) / °° e -'/*P(H, = 0)ds. As in the 
symmetric case P(R S = 0) ~ s" 1 / 2 in d = 1 and s _1 in d = 2, the estimates follow. □ 

Lower bounds. The lowerbounds are through variational formulas. The following is 
the main result of this note and is proved in subsection 2.1. Let e\ and e 2 denote the 
standard basis in M. 2 . 

Proposition 1.3 In d = 2, when Yli^P(^) 7^ an d P = 1/2; there is a constant 
C3 = Cz{p,p) where for all small \, 

/•oo 

2/^(1 -p) / e- xt P(R(t) = 0)dt = ( Vo - p, (A - L)-% - p) p > C 3 log(| log A)|); 
J 

when, more specifically, ^2iip(i) = ce i or ce 2 * s a non-zero multiple of either e\ or e%, 
( Vo - p, (A - L)- 1 ^ - p) p > C 3 \ log A|V2. 



2 Some Preliminaries 

We first give some tools and definitions before going to the proof of Proposition II .31 in 
subsection 2.1. 

Comparison Bound. We compare (/, (A — L)^ 1 ) p with the formula with respect to 
a "nearest-neighbor" operator Lq. Let rrii = • ^ ip{i) for i = 1,2. As the drift of p is 
assumed not to vanish, at least one of the mi's is not zero. Without loss of generality, 
suppose mi 7^ 0. 

Let Lq be the exclusion generator corresponding to nearest-neighbor jump rates 
Po(-) where 

Po(ei) = |mi|,po(e2) = |"i2|> and po(^) = otherwise, when 1712 7^ and 
Po(ei) = |"ii|)Po(±e2) = 1/4, and po(^) = otherwise, when m2 = 0. 

The following is proved in Theorem 2.1 [lOj . 
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Proposition 2.1 There is a constant C4 = C±(p) where 

C 4 - 1 (/,(A-L )" 1 )p < (/, (A - L)~ l ) p < C 4 (/,(A-L )-V 

Duality. Let £ denote the collection of finite subsets of Z 2 , and let £ n denote those 
subsets of cardinality n. Let also $b be the function 

.,<„) - n ^= 



where we take ^0 = 1 by convention. One can check that {^b '■ B G £} is Hilbert 
basis of L 2 (P p ). In particular, any function / G L 2 (P p ) has decomposition 

/-EE k*)** 

n>0 Be£„ 

with coefficient f : £ — > R which in general depends on p. 
Then, for f,geL 2 (P p ), 

<f,0> := </,<7>p = Ef^flW 

and ||f|| 2 := ||/||q = (/,/)«■ Let C n be the subspace generated by local functions of 
degree n, that is functions whose support sets are members of £ n . 

The operators L, S and A have counterparts £, & and 21 which act on "coefficient" 
functions f. These are given in the expressions 

Lf = y,w)( b )*b, sf = E(@f)(£)^, ^ dA f = E^x^- 

Bee Bee Bee 

Let s and a be the symmetric and anti-symmetric parts of p, s(i) = (p(i) +p(—i))/2 
and a(i) = (p(i) - p{-i))/2. Also for B C Z d , denote 

{B\{x}U{y} when x £ B,y B 
B\{y}U{x} when x £ B,y e B 
B otherwise. 

Now, of course, £ = 6 + 21. Moreover, the symmetric part 6 can be computed as 
m(B) = \ E s(y - x)[f(B x , y ) - f(B)}. 

x,yeZ d 

Also, the anti-symmetric part 21 can be decomposed into the sum of three operators 
which preserve, increase, and decrease the degree of the function acted upon: 21 = 
(1 - 2p)2l + 2Vp(1-p)(21+ - 21-). 

(2iof)(#) = E^«(y-^)[f(^)-f(s)] 

y<tB 

(2i+f)(i?) = E-|Ba(j/-^)f(5-M) 

(2l-f)(B) = J2lT B a(y-x)i(BU{x}). 
6 



We note from the expression that in fact 2lo, 2l + and 21 ~ take a degree n function, 
that is say f : (Z rf ) n — > M., into respectively a degree n, n + 1 and re — 1 function. It will 
be helpful to write 21 in terms of its "degree" actions, 

21 = 2l n , n _i + 2l n , n + 21 n,n+l I 

n>0 ^ ' 

where 2l m)n is the part which takes degree m functions to degree n functions. Here, by 
convention 2to,— if = is the zero function. 

At this point, we observe when p = 1/2 that 21 = 2yJ p(l — /o)(2l + — 2l~) as the part 
which preserves degree vanishes here. 

Hi and H_i Spaces. Define, for local functions /, the Hi(P p ) (semi)-norm by 
\\f\\l = (/, {—L)f) p = (/, (-S)f) p . The Hi space then is the completion with respect 
to this norm. With respect to "coefficient" operators, we have the corresponding Hi 
(semi)-norm on functions f supported on £ given by ||f||f = (f, (—£)/) = and 
corresponding completed space Hi. 

Let H-i be the dual of Hi, namely, the completion over local functions with respect 
to norms || • ||_i given by 

\\ff_i = sup{2(f,g) p -\\g\\l} 

9 

= sup{2<f, >-|| 5 ||?} = Hfll^. 
s 

Similarly, we define, for convenience, the notation ||/||^ A = ||f||^ A = (/, (A — S)f) p 
and \\f\\\ x = \\ff_ 1>x = su Pg {2(/, g) p - \\g\\\ x }. 
Then, in this notation, we write for local / that 

</,(A-L)-7>p = sup{2(/,< 7 ) p -|| 5 ||2 A -p 5 || 2 1A } 

9 

= sup{2(f, )-|| S ||2 A -||2l || 2 liA }. 

9 

"Free Particle" Bounds. To analyze these variational formulas, it will be helpful 
computationally to "remove the hard-core exclusion." In other words, we want to get 
equivalent bounds in terms of operators which govern completely independent or "free" 
motions. We follow Bernardin [2]. Let Xn = (Z 2 ) n and note that £ n C Xn- Consider 
n independent random walks with jump rates s on Z 2 . The process xt = {x\, . . . , x") 
evolves on Xn and has generator acting on finitely supported functions 

(6 frcc f)(x) = Y, s(zWx + zej)-ftx)\. 

l<j<n 
z£Z 2 

With respect to finitely supported functions on Xn, let 

(0,^}free = ~f Y] 4>(x)lp(x) 
re! ^— ' 
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be the inner pro duct, and denote the norms ||f||i,free an d ||f||-i,free by 

llf lll.free = h Exe Xn E [fix + ^) " f (x)] 2 

llfll-l.free = sup {2(f,g) free - ||g||? free }. 



Define also the resolvent quantities 



flll,A,frcc 
l-l,A,frcc 



(f, (A — Sf ree )f)f re e 

sup {2(f,0) free - llflllixfreJ. 



Let Q n C Xn be those points whose coordinates are distinct. The following is a part 
of Theorems 3.1 and 3.2 [2] [which simplifies as f = / for / £ C\, and l ie g n f = f for 

f eCn]. 

Proposition 2.2 There exists a constant C§ = C§{n) such that for all functions inC\ 
we have 



a 



- 1 1 1 4:1 1 2 



< 



^5 IN II l,free ^ 11)111 
Also, for all functions in C n (for any n), 

II£||2 ^ r< Il£ll2 



< a 



5||TMl,free- 



|1 M < C 5 \\f\\-l,X. 



free - 



We express now the "free" Hi and H-\ norms in terms of Fourier transforms. Let 
if) be a local function on \n and let ip be its Fourier transform 



-0(si,...,s„) = -= V 



2-Ki{xisi+---x n Sn) 



where s±, . . . , s n E [0, l] 2 . Compute Sf ree from the relation SfreeV' = Sf ree ^ as 



©free^(si) ■■■,Sn) 



where 02(u) = 2 Ylz&z 2 s ( z ) sin 2 (7r(u • z)). 
Hence, 



■0(S1, . . . ,S n ) 



2 

1, A, free 



7 se ([o,i] 2 )™ \ y / 



(si,. . . ,s n )\ 2 ds 



and 



-1, A, free 



,g n )| 2 

6 ([0,1]2)« A + (?2(Sj) 



(is. 
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2.1 Proof of Proposition 11.31 



Let f(rj) = t]q — 1/2. Then, f = (l/2)<5o where 5q is the indicator of the set {0}. 
To prove Proposition 11.31 we find lower bounds on 1 1 <5o 1 1 — i a ■ From Proposition 12.11 we 
will assume L takes form L = Lq. Write now, as p = 1/2, that 



(tfo, (A - = sup(2(5 » 



LA 



The strategy now will be (1) to replace to restrict the supremum on <f> to local degree 
1 functions in C\, and (2) to use the comparison bounds with respect to independent 
particles (Proposition 12.2)) to help bound terms in the variational formula. 

Let now <j> be a degree one function. To simplify notation, let s(±ei) = b\ > 0, 
s(±e 2 ) = 62 > and a(e%) = — a(— e\) = a\ ^ 0, a(e2) = — a(— e 2 ) = 02- Note now 
that 210 takes form 210 = 211,20. More specifically, (2li )2 0)(x, y) = a(y — x){4>{x) — 4>{y)) 
is supported on distinct two-tuples and can be written as 

ai((p(x) — 4>(x + ei)) when y = x + e± 
(2li,20)( x >2/) = { a 2(4>(x) - 4>(x + e 2 )) when y = x + e 2 

otherwise. 



Using Proposition 12.21 we have for some constant Cg that 



l,A,free|- 



(6o, (A -£)"%> > C7 6 sup ^2<(yo J 0>-||0||l A , frBe -||Sli I 20||ii,A,free^ (2.2) 

local in Ci 

Now, it is a calculation to find that 

= ^ E e 2 ™^+^(2li, 2 0)(*,y) 
z 

= — - ^=4>(s + t)[2a\ sin(27rsi) + 2a 2 sin(27TS2) 
a/2 

+2ai sin(27rti) + 2a2 sin(27ri2)]- 
Then, the expression in brackets in (|2.2|) becomes in Fourier terms 

s)-(\ + e 2 ( s )M 2 )ds 

[0,1] 2 ' 

1 f ( EL ^ sin(27r Sj ) + 2a, sin(2^)) 2 ^ + t)? 



2 i ([ o,i]2)2 x + e 2 (s) + e 2 (t) 

We now change coordinates in the second integral: 

,u + v w + z u — v w — z 



{si,s 2 ,h,t 2 ) = (- 
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(whose Jacobian determinant in absolute value is 1/4). The region [0, l] 4 is mapped 
to D 2 where D is a planar diamond with vertices (0, 0), (1, —1), (1, 1), (2, 0). Let 

2 . U + V o, W + Z 

j{u,v,w,z) = 2oi sin (tt — - — ) + 2o 2 sin (ir — - — ) (2-3) 

+2&i sin 2 (7r— ^— ) + 26 2 sin 2 (vr— ^— ) 

= 4&i sin 2 (vr(u/2)) cos 2 (vr(?;/2)) + 46i sin 2 (vr(?j/2)) cos 2 (vr(u/2)) 

+46 2 sin 2 (vr(u;/2)) cos 2 (vr(z/2)) + 4& 2 sin 2 (7r(z/2)) cos 2 (vr(u;/2)) 

and 

v(u, v, w, z) = 4a 2 sin 2 (7ru) cos 2 (7re) + 4a 2 . sm 2 (irw) cos 2 (ttz) 
+8aia 2 sin(7ra) cos(7ru) sin(7ui;) cos(7rz) 
< 8af sin 2 (7ru) cos 2 (tt?;) + 8a 2 . sin 2 (7ra) cos 2 (7rz) . 

The second integral is rewritten as 

v(u, v, w, z) 




D JD 



X + ^(u,v,w,z) 



u,w)\ dudvdwdz 



, , 8a? sin 2 (-7ra) cos 2 (7rw) + 8an sin 2 (7rw) cos 2 (7rz) -~ 2 

< / / : ; : \<p(u,w)\ dudvdwdz. 

~ 1 1 X + j{u,v,w,z) 1 y n 




D JD 



By changing variables and adding and identifying some parts of the region of inte- 
gration, it is not difficult to see that the last integral reduces to 



4 / / / / — — — - — '- — — i — -\(j)(u,w)\ z dvdzdudw 

Jo Jo Jo Jo X + j{u,v,w,z) 

= 4/ / [8ajsm 2 (7ru)Fl(u,w) + 8alsm 2 (7rw)F 2 (u,w)]\^(u,w)\ 2 dudw 
Jo Jo 

where 

r,\t \ f 1 f 1 cos 2 (irv) , , , . f 1 f 1 cos 2 (irz) , , 

l< x (u,w)= / / -dvdz and F x (u,w) = / / 7 -dvdz. 

Jo Jo X + j(u,v,w,z) ' J J X + j(u,v,w,z) 

Let now C7 = 2(6f + b 2 ,) + 16(a 2 + a 2 )- For general rates, substituting into (|2.2[) . 
we obtain C f r 1 (5o ) (A — ii)^ 1 ^) greater than 

sup^ / / 0(«,u>) 
</> Uo Jo 

-(A+ C 7 (sin 2 (7ru) + sin 2 (7ro))(l +F^(u,w) + F^(u,w)))\(j)(u,w)\ 2 dudw \ 
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where the supremum is on <j) local, or without loss of generality on L 2 {£{). When 
a 2 = 0, we have the lower bound 

sup( / / $(u,w) (2.4) 



— (A + CV(sin 2 (7ru) + sin 2 (7r?x;)) + Cj sin 2 (iru)Fl(u, w))\(f>(u, w)\ 2 dudw^j . 

We now concentrate on the general rates case. By optimizing on <j) we get the lower 
bound 

1 f 1 f 1 dudw 

4 J J A + C 7 (sin 2 (vra) + sin 2 (vru;))(l + F^u,w) + F 2 (u,w)) 
with optimizer 

d>(u,w) 



2 A + C 7 (sin 2 (vra) + sin 2 (vrw))(l + F^(u, w) + F 2 (u, w)) 

which is the transform of a real function as (f>(u, w) = 0*(1 — u, 1 — w) (note F^(u, w) = 
FUl — u,l — w) for i = 1,2 by observing Q2.3|) and changing variables v — ► 1 — v and 
z->l-z). 

We now bound Fl(u,w) + F 2 (u,w) for \u\, \ w\ < 1/2. Since cos(x) is decreasing for 
< x < tt/2 and sin(aj) > (2/ir)x for < x < tt/2, we have 

/ ^9 / \ f 2dvdz 

Fl{u,w) + Fl{u,w) < / 



[([i/2,i]x[o,i])u([o,i]x[i/2,i]) A + 2biv 2 + 2b 2 z 2 

V2 fl/2 2dvdz 

+ ' 



Jo 



X + 2b 1 (u 2 + v 2 ) + 2b 2 (w 2 + z 2 ) 



< C 8 + 



n / 2 f 1 2sdsda 



A + b(u 2 + w 2 + s 2 ) 



vr ^ r A + 6(n 2 + ^ 2 ) + 5 1 
= 6 8 + -= log - — 

< C 9 + ^ b \log(X + b(u 2 + w 2 ))\. 

where b = 2min{6i, b 2 } and Cg = C%{b\, b 2 ), Cg = Cg(fei, b 2 ) are constants. 
Hence, as sin(x) < x, we can bound (|2.5[) below by 

1 r 1 / 2 r 1 / 2 

4i 7o A + 4C 7 ^ 2 (n 2 + «; 2 )(l + C 9 + ^|log(A + 6(« 2 + u; 2 ))|)' 

We have with respect to constants C±o, C\\ that 4 times (j2.6[) is greater than 

1/2 rdr rl/(2VX) rdr 



(2.6) 



A + Cior 2 (l + |log(A + dor 2 )|) J 1 + C 10 r 2 (l + | log A(l + C 10 r 2 )\) 
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1 /• 1 /(2v / A) dr 

> 



C\\ J i r \ log Ar 2 | 

1 r 1/2 dr 



Cll iv^ r|logr 2 |' 

This last expression is order | log(log A)|. To get the larger expected order of | log A| 2//3 , 
it seems one would need to optimize also over higher degree functions in (|2.2|) . 
We note in the case 02 = 0, we bound (|2.4j) by 

1 f 1 / 2 f 1 / 2 dudw 

i/o Jo A + 4<^vr 2 (u 2 + w 2 ) + C' 7 u 2 {C^ + f b \ log (A + b(u 2 + w 2 ))\) ' 

Following closely the sequence to bound the second-class particle variance in d = 2 (cf. 
p. 470 0), we observe | log(X +u 2 +w 2 )\ < \ log(A+w 2 )| for A small and < u, w < 1/2. 
And so, we obtain a lower bound on order 

V 2 rV 2 dudw 



A + u 2 + w 2 + u 2 \ log (A + w 2 )| ' 
With substitution u = y(l + | log (A + w 2 )\) the above expression is bounded below by 

- 1 / 2 TV 2 rfy^ 



\ + y 2 + w 2 



(l + |log(A + u; 2 )|)- 1 / 2 . 



Changing to polar coordinates and restricting tt/6 < a < vr/4, we get a lower bound 
on order as in 

1/20 rdr 

T^|log(A + r 2 )|- 1/2 > C 12 |logA|V 2 
A + r z 

for a constant Cyi- □ 
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